On the Instability of the Lee- Wick Bounce 
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It was recently realized [l[ that a model constructed from a Lee- Wick type scalar field theory yields, 
at the level of homogeneous and isotropic background cosmology, a bouncing cosmology. However, 
bouncing cosmologies induced by pressure-less matter are in general unstable to the addition of 
relativistic matter (i.e. radiation). Here we study the possibility of obtaining a bouncing cosmology 
if we add radiation coupled to the Lee- Wick scalar field. This coupling in principle would allow the 
energy to flow from radiation to matter, thus providing a drain for the radiation energy. However, 
we find that it takes an extremely unlikely fine tuning of the initial phases of the field configurations 
for a sufficient amount of radiative energy to flow into matter. For general initial conditions, the 
evolution leads to a singularity rather than a smooth bounce. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

Both Standard and Inflationary Cosmology Q suf- 
fer from the initial singularity problem and hence cannot 
yield complete descriptions of the very early universe. If 
one were able to construct a non-singular bouncing cos- 
mology, this problem would obviously disappear. How- 
ever, in order to have a chance to obtain such a non- 
singular cosmology, one must either go beyond Einstein 
gravity as a theory of space-time (see e.g. [J| for an early 
construction) , or else one must make use of matter which 
violates the "null energy condition" (see Q for a review 
of both types of approaches). 

Interest in non-singular bouncing cosmologies has in- 
creased with the realization that they can lead to alterna- 
tives to inflationary cosmology as a theory for the origin 
of structure in the universe. A specific scenario which can 
arise at the level of homogeneous and isotropic cosmol- 
ogy is the "matter bounce" paradigm which is based on 
the realization [(| @] that vacuum fluctuations which exit 
the Hubble radius during a matter-dominated contract- 
ing phase evolve into a scale-invariant spectrum of cur- 
vature perturbations on super-Hubble scales before the 
bounce. The key point is that the curvature fluctuation 
variable £ grows on super-Hubble scales in a contract- 
ing phase, whereas it is constant on these large scales in 
an expanding phase. Since long wavelength modes exit 
the Hubble radius earlier than short wavelength ones, 
they grow for a longer period of time. This provides a 
mechanism for reddening the initial vacuum spectrum. 
It turns out that a matter dominated contracting phase 
provides the specific boost in the power of long wave- 
length modes which is required in order to transform a 
vacuum spectrum into a scale-invariant one. Studies in 
the case of various non-singular bounce models Q have 
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shown that on wavelengths long compared to the dura- 
tion of the bounce phase, the spectrum of fluctuations is 
virtually unchanged during the bounce. Thus, a scale- 
invariant spectrum of curvature fluctuations survives on 
super-Hubble scales at late times. 

Provided that the bounce can occur at energy scales 
much below the Planck scale, non-singular cosmologies 
solve a key conceptual problem from which inflationary 
cosmology suffers, namely the "Trans-Planckian" prob- 
lem for fluctuations [1, [l(J : If the period of inflationary 
expansion of space lasts for more than 70-ff -1 , where H 
is the Hubble expansion rate during inflation (in order 
to solve the key cosmological mysteries it was designed 
to explain, inflation has to last at least 50i/ _1 ), then 
the physical wavelengths of even the largest-scale fluc- 
tuation modes we see today will be even smaller than 
the Planck length at the beginning of inflation and thus 
in the "zone of ignorance" where the physics on which 
inflation and the theory of cosmological perturbations 
are based, namely Einstein gravity coupled to semiclas- 
sical field theory matter, will break down. In contrast, 
in a non-singular bouncing cosmology the wavelength of 
modes which are currently probed by cosmological ob- 
servations is never much smaller than 1mm (the phys- 
ical wavelength of the mode which corresponds to our 
current Hubble radius evaluated when the temperature 
of the universe was 10 16 GeV) and hence many orders 
of magnitude larger than the Planck length. Thus, the 
fluctuations never enter the "trans-Planckian zone of ig- 
norance" of sub-Planck-length wavelengths. 

Possibly the simplest realization of the matter bounce 
scenario is the "quintom bounce" model fill ) and is ob- 
tained by considering the matter sector to contain two 
scalar fields, one of them (the "ghost field") having the 
"wrong" sign of the kinetic action. The potential of the 
ghost scalar field also has the opposite sign to that of 
regular scalar fields such that in the absence of interac- 
tions, the ghost field has a classically stable minimum. 
As has been noticed in [l[ , such a quintom bounce model 
also arises from the scalar field sector of the "Lee- Wick" 
(LW) Lagrangian [ijj which contains higher derivatives 
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terms. 

The quintom and Lee- Wick bouncing cosmologies are 
obtained in the following way [l|, [ll|: We begin in the 
contracting phase with both the regular and the ghost 
scalar field oscillating homogeneously in space about 
their respective vacua. We assume that the energy den- 
sity is dominated by the regular matter field, and that 
hence the total energy density is positive. Once the 
amplitude of the regular scalar field exceeds the Planck 
scale, the field oscillations will freeze out and a slow- 
climb phase will begin during which the energy density 
of the field only grows slowly (this is the time reverse of 
the slow-roll phase in scalar field-driven inflation). How- 
ever, the ghost field continues to oscillate and its en- 
ergy density (which is negative) continues to grow in 
absolute value. Hence, the total energy density drops 
to zero, at which point the bounce occurs, as has been 
studied both analytically and numerically in the above- 
mentioned works. Note that the energy density in this 
bounce model scales as matter until the regular scalar 
field freezes out. 

A major problem of bouncing cosmologies realized with 
matter which scales as a~ 3 as a function of the scale fac- 
tor a(t) is the potential instability of the homogeneous 
and isotropic background against the effects of radiation 
(which scales as a -4 and anisotropic stress which scales 
as a~ 6 16!]. If we simply add a non-interacting radiation 
component to the two scalar field system, then unless 
the initial energy density in radiation is tuned to be ex- 
tremely small, then the radiation component will become 
dominant long before the bounce can arise, and will pre- 
vent the energy density in the ghost field from ever being 
able to become important, resulting in a Big Crunch sin- 
gularity. Similarly, unless the initial energy density in 
anisotropic stress is very small, it will come to dominate 
the energy density of the universe long before the bounce 
is expected. The anisotropics will destabilize the homoge- 
neous background cosmology, and will prevent a bounce. 
Note that at the quantum level, there is an additional 
severe problem for bounce models obtained with matter 
fields with ghost-like kinetic terms, namely the quantum 
instability of the vacuum (see e.g. [3). 

In this paper we will focus on the radiation instability 
problem. For the purpose of this discussion we will sim- 
ply assume that anisotropic stress is absent. In a recent 
paper, two of us studied the possibility that a bounce 
could arise if radiation is supplemented with Lee- Wick 
radiation jTEj . However, we showed that this hope is not 
realized: the addition of Lcc-Wick radiation does not pre- 
vent the Big Crunch singularity from occurring. In the 
presence of radiation, the only hope to obtain a bounce 
is to introduce a coupling between radiation and ghost 
scalar field matter which could effectively drain energy 
density from the radiation field and prevent the energy 
density of radiation from becoming dominant. Here we 
study this possibility. However, at least for the specific 
Lagrangian which we consider, we find that a bounce 
only emerges for highly fine-tuned phases of the fields 



and their velocities in the initial conditions. 

The paper is organized as follows: In Section II, we in- 
troduce the model we study, namely the scalar field sector 
of Lee- Wick theory coupled to radiation, and write down 
the general equations of motion. In Section III we set 
up the equations of motion linearized about the bounce 
background, treating the entire radiation field as an inho- 
mogencous fluctuation. In particular, we study the dif- 
ferent terms which contribute to the energy-momentum 
tensor and identify those which could assist in obtaining 
a non-singular bounce. In Section IV we study the solu- 
tions of the perturbed equations of motion, and in Sec- 
tion V we analyze the evolution of the different terms in 
the energy-momentum tensor, identifying the conditions 
which would be required in order to obtain a non-singular 
bounce. We have also evolved the general equations of 
motion for the two inhomogeneous scalar field configu- 
rations and the classical inhomogeneous radiation field 
in the homogeneous background cosmology. Section VI 
summarizes some of the numerical results. Both the ana- 
lytical and numerical results confirm that we need unnat- 
ural fine-tuning of the initial conditions in order to obtain 
a non-singular bounce. In the final section we offer some 
conclusions and discussion. 



II. THE MODEL 

The Lee- Wick scalar field model coupled to electro- 
magnetic radiation is given by the following Lagrangian: 

c = -I^a^ + ^(aW-^ 2 ^-y(0) 

--F^ -j^^F^Fn , (1) 

where m is the mass of the scalar field cf>, and V(4>) is its 
potential. Here we adopt the convention that 

ds 2 = -dt 2 + a 2 {t)(dx 2 + dy 2 + dz 2 ) , (2) 

where a(t) is the scale factor of the universe. Since it is a 
higher derivative Lagrangian in 0, the scalar field sector 
contains an extra degree of freedom with the "wrong" 
sign kinetic term and with a mass set by the scale M. 
We choose m < M < mpu where mpi is the Planck 
mass, since we want the regular scalar field to dominate 
at low energies, but at the same time we do not want 
to worry about quantum gravity effects. The second line 
of the Lagrangian (JTJ) contains the kinetic term of the 
radiation as well as the coupling term, where we assumed 
both for the sake of generality and because of foresight 
that the radiation field couples not only to the scalar 
field <j> itself, but also to the higher derivative term. The 
electromagnetic tensor, F^, is related to the radiation 
field through the usual definition 

= VnA v - V„A M , (3) 

where V M is the covariant derivative. 
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It is convenient to extract the extra degree of freedom 
as a separate scalar field. To do this, we use the field 
redefinitions 

<t> = 0i - <h , 

fa = d 2 fa/M 2 . (4) 
The Lagrangian ([1} then takes on a simpler form: 

C = -^fad^fa + \d„fad»fa - \rr?<& + ^M^l 

-^F^-fWufaF^Fn , (5) 

where we have chosen the potential to be zero. In this 
new form, the Lagrangian describes two massive scalar 
fields with one of them (i.e., fa) behaving like a"ghost", 
and both of them coupled to the radiation field. 

The coupling term f(fa, fa, F^ U F^ V ) should in prin- 
ciple be arbitrary, however, in this paper we will take a 
specific form for convenience. The form will be: 

V") = \{c4>l + d4>l)F^F^ , (6) 

where c and d arc coupling constants which have 
mass dimension —2. The interaction terms are non- 
renormalizable. To make sure that such terms could 
be thought of as arising from an effective field theory 
which is consistent at the bounce, we must make sure 



that the coefficients are chosen such that the contribu- 
tion of the interaction term to the Lagrangian density is 
smaller than that of the other terms. This must be true 
even at energy densities at which the bounce occurs in 
the pure scalar field model. It is easy to see that this 
condition will be satisfied if the coefficients c and d are 
both of the order to~ ; 2 . 

It is the purpose of this paper to study the effects which 
these coupling terms have on the dynamics of the sys- 
tem. We know that in the absence of coupling, i.e. when 
c = d = 0, a bounce will only occur if the initial radiation 
energy density is tuned to a very small value compared to 
the scalar field energy density. This is because the posi- 
tive definite energy density of radiation will scale as a -4 
which is faster than that of the scalar fields, in particular 
the ghost scalar field. Generically, it will dominate the 
energy of the universe after some amount of contraction, 
it will prevent the ghost scalar field energy density from 
catching up and will thus prevent a bounce, leading to a 
Big Crunch singularity instead. With non-vanishing val- 
ues of c and d, however, the scalars are in principle able 
to drain energy from the radiation. 

From the Lagrangian ([5]), one can obtain the stress- 
energy tensor by varying the action with respect to 
the metric g^ v . In the hydrodynamical limit, we can take 
Tp, v to be of the form of diag{p, a 2 (t)p 1 ,a 2 {t)p2,a 2 {t)pz\ 
where p and p are energy density and pressure, respec- 
tively. For now, we consider the general form which is: 



Tp V = g^C + d^fadyfa - d^fadufa + (1 - c<j>\ - dfa^F^xF* 



9vu [ y - ^-pdifadifa 
+d^fad v fa - d^fadvfa + (1 - aj>\ 



2 ±2 
-TO fa 



- dfa^F^Fl 



2 d r 



2+2 



1* 



c</> 2 - dfa 2 )F 2 



(7) 



Since we will be studying the contribution of plane wave 
perturbations of the scalar fields and we will treat ra- 
diation as a superposition of waves, we kept the space- 
derivative terms. 

By varying the Lagrangian with respect to the matter 
fields fa, fa and A^, we also get the equations of motion 
for all three fields: 

Ufa - (to 2 - C -F 2 )fa = , (8) 

Ufa - (M 2 + ^F 2 )fa = , (9) 

(1 - c4 2 - dfa 2 )(d u F ilv + 3HF»°) 
-2{cfad y fa + dfad y fa)F^ = (10) 

which will be analyzed in detail in the rest of the paper. 



III. DYNAMICS 

Since the equations of motion are nonlinear, we can- 
not work in Fourier space, and use plane wave solutions. 
However, we are interested in how initially small amounts 
of radiation build up and possibly transfer their energy 
to scalar field fluctuations. We treat radiation as a super- 
position of fluctuations. Therefore it makes sense to lin- 
earize our equations about the homogeneous scalar field 
background. Thus, we make the following ansatz for the 
scalar fields: 

fa(t,z) = 4° ) (t)+6^ 1) (t,z)+6 2 ^ 2) (t) (11) 

fa (t, z) = 4 )] (t) + e4 1] (t,z) + 6 2 4 2) (t) , (12) 

where the expansion parameter e is taken to be much 
smaller than 1 [l?} • The first term on the right hand side 
of each line, i.e. fiili't) correspond to the background 
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fields, the terms <j>\ 2^t x ) are the fluctuations, and the 

second order terms 4>\ ' 2 (t) describe the back-reaction of 
the fluctuations on the background and can be computed 
from the leading second order corrections (averaged over 
space) of the equations of motion [l8| . 

To simplify the analysis, we describe radiation in terms 
of plane waves in a fixed direction (which we take to be 
the z direction). Without loss of generality we can re- 
strict attention to one polarization mode which we take 
to be the electric field in the x direction and the mag- 
netic field in the y direction. In this case, the only non- 
zero components of the field strength tensor are F 01 and 
F 13 . Using the temporal gauge where A = 0, we find 
that only the first component of the gauge field is non- 
zero. For a single wavelength fluctuation we can make 
the ansatz 

Ai(k,t) = f{t)cos{kz) = 1 {k,t) , (13) 



or, equivalently, 

A 1 (k,t) = a(t)- 2 ~/(k,t) . (14) 

Since in the linearized equations of motion the Fourier 
modes are independent, we can consider <fr^ and <j)^ 
also to be plane waves propagating in z direction, so they 
depend only on z and t. 

With Eqs. (|1HI14[) in hand, we can write down the 
energy densities of the various fields at each order in per- 
turbation theory. 

A. The stress-energy tensor 

First of all, we insert the above perturbative ansatz 
for the fields into the stress-energy tensor of the system. 
From the general expression ([7]) for T^ v we get: 



1 



d z (bid z 4>i 



1 



2/2 ^2 

2 m ^ ~ T 



2a 2 

- d^fcdvfc + (1 - af>l - 



2a 2 



-M 



1 



(l - c4>l 



\)F 2 



(15) 



The 00 component of Eq. (|T5|) denotes the energy den- 
sity of the system: 

p = \^\ + %^ + m 2 ^)-\^ 2 + %4>l + M 2 ^ 2 ) 

2 a z 2 a z 



+ (1 - ccjij - dq> 2 )(^ + F 0X F X ) , 
so at each level in perturbation theory we have: 



(16) 



We can similarly obtain the pressure of the system from 
the ii components of Eq. (fl~5"j) . Note that due to the 
anisotropy in T^ v caused by the gauge field as well as 
by the anisotropic fluctuations of the scalar fields, the 
pressures in the three directions are no longer identical. 
The pressure in each direction can be written as: 



,(2) 



p<°> = l^ + rn^-^' + M*^) (17) 



1 



+m 2 4 0) 4 1) ) , 



-.a 



(18) 



m (f>\ 



\l-%<t,l-m^\) 
a z 



-(1 



F 2 



4< 

FxF, x , 



a 2 



J2 - -0 - M 



(20) 



+M 2 4 1)2 + M 2 4°V 2 2) ) 

+ (l-c0l O)2 -d^ 2 )(^- + F 0X F x ) , 



(i) . ;(0) ;(2) 



(i) 



2^2 



1 



(1) 



4< O) 0< 2 > 



2^(1)^ 



m <fi 



with no summation over the index i. From this formula, 
we can see that at both zero-th and first order, the pres- 
sure is isotropic: 



+m 2 ^°M 2) ) - ^ +02 O V2 2) + ^F 



+M 2 ^\M 2 ^Ui 2) ) 



1 - a 



2 , , k 



7 



'2a 



2a 2 ' 



(0) 



j(0) 



™2 ,(0)" 



2 ), 



4°> 



M 2 ^ 0) )(21) 



WW 



-M 



2^(0U(lh 



(22) 
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while the second order pressure for each direction gives 
1 , ,m 2 ;«r> U2) k 2 



P 



(2) 



3^1 



2j,W 



m 0{ 



w^vi 2) )-i(4 1)2 +4 0) 4 2) - 



-A/ 



,(o) ,(2), . di<f>idi<f>i d l (j) 2 d t (j)2 



a- 



-(l-c^) 2 -^ 2 )^-^), 



(23) 



where i 


= 1,2,3. 






(2) 

We can thus obtain every component of p\ . 




( 2 ) 


1-2 ;„2 n 


2 i(l) 2 
TO 0^ 






+m 2 ^ O Vi 2) )-^(4 1)2 + 2 O) « 




A (l) 2 
:; 2 












+ (l-<) 2 -M° )2 )(^7 2 - 


2fl2 j ' 


(24) 


( 21 




2 (l) 2 
TO 0} 






1 • 2 

2 /Col ,f21\ 1 / ;(o) 

+m 2 <Vr) - 2(4 J +4 « 


1(2) fc2 


,(i) 2 




i i. f 2j,(l) 2 , »,r2 .(0) ,(2)n 








-(l-<) 2 -M 0)2 )(^7 2 - 


2a2 j ' 


(25) 


(2) 
2>3 = 




TO 2 !^^ 








or 


4 1)2 




+M 2 4 1)2 + M 2 4 0) ^ 2) ) 








+ (l-<) 2 -M° )2 )(^7 2 - 


7 2 


(26) 


and the 


average is: 






„(2) _ 


(2) , (2) . (2) 
Pi +P2 +P3 

3 




(27) 




= i(^) 2 +^)^)-^) 2 + 


TO 0^ 






+m 2 ^^ 2) )-i(4 1)2 + # 


1(2) ^ 
02 





+M 2 2 1} +71/ 



2 4 0) 4 2) ) 

,(o)V fc2 ..2 , 7 2 



+(l-c0r -d^ )(^T7 2 + ^). (28) 

From the above, we can also see that in order to an- 
alyze the behavior of the energy density up to second 
order, we need to know the evolution of scalar fields up 
to second order as well as that of the gauge field up to 
first order, while the behavior of the gauge field to second 
order is not required. 



It will be useful in the following to separate the contri- 
butions to the energy density and pressure in a different 
way, namely, 

i) the contribution from the background homogeneous 
part of the scalar fields, 



(0) 



n h - n (0) 
P4> —Pcf, , 



(29) 



ii) that of the scalar field perturbations (in slight abuse 
of notation we call this the "inhomogcneous" term) , 



iii) the contribution of the gauge field, 
1 h 2 

/ ■ 2 2 \ 

p ° = ii (72 + ^ 2) ' 

and iv) the contribution of the coupling term, 

^2 m\2 

Pc = 



8a 2 



(30) 

(31) 
(32) 

<5> Pg , (33) 



Pc 



where in the last equation we define $ to be the quadratic 
combination of the two fields: 



$ = (c^ ] + d^" )/2 



(0)' 



(35) 



From the above we can deduce the equation of state 
parameter for each part: 



, (36) 



inh 



< )2 - 


to 2 ^ 2 - 


4 0)2 


+ Af 2 4°) 2 


4 0)2 + 


2^.(°) 2 


4 0)2 


-m 2 4 0)2 


(p [i) 










p(°) 2 , 


) 





/gW p(0)p(2) 

^ 9 (0) 2 o (0) 2 I ' 



1 

3 



P 1 



(38) 



From the equations above, we see that for positive val- 
ues of the constants c and d, the coupling of the scalar 
field with the gauge field will give rise to a contribution 
p c to the energy density which has the same equation of 
state but opposite sign to that of the gauge field. There- 
fore the coupling can help drain energy from the gauge 
field. It is because of this mechanism that we might hope 
to achieve a cosmological bounce in the presence of ra- 
diation. A first indication on whether a bounce might 
occur can be obtained by considering the scaling of each 
contribution to the energy density as a function of the 
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scale factor a(t). To find these scalings, we need the 
time dependence of the linear and quadratic contribu- 
tions to each field. Therefore, we need to solve the matter 
field equations of motion. In the following subsection, we 
present the equations for the fields at each order, while 
the solutions and detailed analysis will be performed in 
the next sections. 



B. Equations of Motion 



is used rather than the cosmic time. Additionally, to ex- 
tract the dependence on the scale factor, it is convenient 
to use the following two variables: 



ui(rf) = a{>q)4>i{rj), u 2 {r]) = a{rf)4>i{ri) 
.«/„• 



(43) 



Hereafter, we will use uy (i = 0,1,2, j = 1,2) to de- 
note the i-th order perturbation of the j-th scalar field. 
Moreover, for simplicity but without loosing generality, 
we can parameterize the scale factor a(t) as 



Keeping in mind the ansatze for A^, fa and fa , their 
equations of motion at each order can be obtained from 
©, © and OH): 
a) At zero-th order: 



1(0) 



■ 3ff^ 0) + m 2 <^ 0) = , 



3H^ 



M 



(39) 



= 



Note that there is no equation at this order for be- 
cause it is of first order in e. 
b) At first order: 



^ + + (£ + m 2 )^ 1 ' 



M 





: 



(40) 



z0\ o)2 - dfaf )2 )(d„F^ + 3HF p0 ) 



(1-00 



(o) a ^(o). 



= 



Making use of Eqs. (|T3")) and (|T4"|) . the equation for the 
gauge field can also be rewritten as: 



-2(c0r^ o) +^ o) 4 o) )7=o. a 



k 2 



c) At second order: 



mrq>\ ' 



f < F fitll F tlv > Vl 



(0) 



M 



2^,(2) 



4 < F UU F^ > 



(41) 



(42) 







Here, pointed parentheses indicate spatial averaging 
(since we are only focusing on the zero mode of the second 
order field fluctuations). We also neglected the effect of 
Hubble friction since it does not give an important con- 
tribution for the second order fluctuations. 



IV. THE GENERAL SOLUTION 

In this section we will solve the equations of motion 
([39| . (|40| and (|42|) to see if and how a bounce will hap- 
pen. 

It is usually useful to perform the analysis in the con- 
formal frame where the conformal time i] = Ja _1 (i)di 



i(rj) = a Q t p = aoM l -» , 



with 



3(1 + w) 



(44) 



(45) 



where ao and w are the initial value of the scale factor 
and the equation of state of the universe, respectively. 
This is a self-consistent assumption when w is nearly a 
constant. The evolution of w in our case will be shown 
numerically in Section VI. 



A. Solutions for <jy± ' and cf)^ 



Using the parametrization (|44p . the equations of mo- 
tion at zeroth order of the two scalar fields become: 



(46) 



where a prime denotes the derivative with respect to con- 
formal time rj. Their solutions are: 



,(o) 



,(()) 



#±i=a((l-p)amH), (47) 



\V\H ±1 ^((1 - p)aM\ V \) , (48) 



where H ± is P represents the (± 2 p )-th order Hankcl 
function. Far away or close to the bounce, i.e. for 
a\rj\ 3> m _1 ,M _1 and a\r]\ <C m _1 ,M _1 , respectively, 
the approximate solutions are: 

1) Oscillations for large values of the scale factor a\rj\ 3> 
m _1 ,M _1 : 



cos((l - p)a m\r)\ l -v 

y (1 — p)ira m v 

+0{ 0) ) , (49) 

((1 -p)a M\ri\^ 



(1 — p)ira M 



cos 



+4 0) ) , 



(50) 



(where 9\ and 62 are phases set by the initial conditions). 
In terms of the non-rescaled fields fa one obtains damped 



7 



(or anti-dampcd) oscillations (depending on whether we 
are in an expanding or a contracting period) 



^ 0) ~ M^ry 



i 0) ) 



(1 — p)Tra§m 



>s((l - p)a m\ri\ 



(51) 



(1 -p)ira£M 



cos((l — p)a M\r]\ 1 -p 



(52) 



2) "Frozen" evolution for small values of the scale factor 
o|ry| < m~ l , M -1 : 

/ \ 1-3p l-2p 

(0) ((1 -p)a m) 2 |t?| 



r( Ml f p)) 

{(l-p)a m) 2 |?7| i-p 

(0) ((i-p) fl0 M)^|^ 



(53) 



r(toi) 
r(ifH) 



(54) 



from which it follows that the non-rescaled fields 
evolve as 

,(0) ({l-p)a m) 2 \rj\ i-p 



a r(M) 

3p-l 

((l-p)a m) 2 

+ « r(i^) ' 

,(o) ((l-rtapM) 1 ^^!^ 

2 ~ a r(^) 



(55) 



((l-p)a Af)^ 



(56) 



(.(0) 



from which we can see that the last term of cp\~' is a 
constant mode while the first term is a varying one. De- 
pending on the value of p (or equivalently w) the varying 
mode could be growing (for p > 1/3 or —1 < w < 1), in 
which case it becomes dominant, or decaying (for p < 1/3 
or w > 1/w < —1), in which case it becomes subdom- 
inant. We can usually neglect the decaying part of the 
fields. 



B. Solutions for <j>Y and tf>^ 



Following the steps performed in the last subsection, 
we can also get the solutions for the first order compo- 
nents of the scalar fields. Using the equations (|40|) for the 



first order perturbations we obtain the following equa- 
tions of motion for ujj 1 and «2 : 



i 1) " + (^ + a^^-^)4 1) =0. 



(57) 



Depending on the value of k, we obtain different ap- 
proximation solutions. For wavenumbers large compared 
both to the Hubble radius and to the mass term, we ob- 
tain oscillatory solutions with fixed amplitude. 

Considering now modes which arc still sub-Hubble (i.e. 
k\r]\ > 1) but for which the mass term dominates over the 
contribution of the field tension (i.e. the term involving 
k), we can neglect both the k 2 term and the term involv- 
ing ri~^p^3 ■ The simplified equation for these modes 



(l)" „ „ 2p (1 ) 



(58) 



+ a^M z r]T^u\ ' = , 



whose solutions are: 



v\Hi^p((1 -p)a m\r]\) 



(1 — p)irmao 



cos 



((1 -p)a m\i]\ 



(59) 



v\H^((l-p)a M\r)\) 



(1 — p)irMao 



i(l - p)a M\r ] \ 



(60) 



For modes outside the Hubble radius (kr/ <C 1), we 
have: 



uP" + (ogmV?* - ^)u? = 



(61) 



which have the same form as Eq. (|46[) so their solution 
will be the same as given in Eqs. (|53[) and (|54"|) . 

We have thus seen that the first order solutions for the 
scalar fields scale the same way with \rj\ as the zeroth or- 
der solution. This is because in the small \r/\ region where 
the a"(t)/a(t) term dominates over the other ones, the 
equations for first order and zero-th order modes are al- 
most the same. Thus, unless the energy density in the 
uj modes dominates at the initial time, it will never 

dominate over the background contribution from the 
terms. Thus we can conclude that the first order fluctu- 
ations of scalar fields will not prevent the bounce. 
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C. Solution for the gauge field 7 

In this section, we will analyze the gauge field 7 which 
is also considered to be of first order. The equation (|4Tj) 
can directly be transformed to conformal frame as: 



2(c6 {0) cb {0) ' + dcb {0) 6 {0) ') 



(oy 



(62) 



Since the coefficients c and d are small, we can take the 
last term to be a source term. In a first order Born 
approximation, we can write the total solution as 



7 ~ 70 + Sj , 



(63) 



where 70 is the solution for the homogeneous equation 
obtained by setting c = d = 0, while S"f is the leading 
correction term obtained by inserting 70 into the source 
term (the last term in (f6"2"]) ). 

The zero-th order (homogeneous) equation is easily 
solved and gives 

7 O ~cos(fc|?7|+0 7 ) . (64) 
For the first order equation, it is convenient to define 

2(c6 {0) 6 { ° y + d6 {0) 6 {0) ') 
p() __^L+*i2j (65) 

i-c4 0)2 -d4 0)2 



so that the equation becomes 

£7" + k 2 5-y + P{ri)i = , 



(66) 



where we neglected the small term P(rj)8j'. Inserting 
the solution of 70 we get the following equation for 
^7: 



fry" + k 2 8 1 = -P{vH 



(67) 



We are interested in the scaling of 8^ as a function of 
time. For this purpose, we need to work out the scaling 
in time of the source term in (|67[) . Since the solutions 
for </>^ and (f> 2 scale differently in time in the two time 
intervals discussed in Subsection pV A[) . it is necessary 
to analyze these two intervals separately. 

For times obeying a\r)\ ^> m , M , then by differen- 
tiating (f5"Tj) and (f52"j) with respect to r\ we have: 



3p 



(1 — p)ira,Q m 



(68) 



p-1 



s((l — p)aom\i]\ l -p 



,(0h 



+aom\r)\ i-p sin((l — p)a m\r)\ 1 ~p +6^' )] , 

(69) 



1 1 p + 2 
\rj\ 2 <p-!> 



(1 -p)na$M 



: f — ^- cos((l - p)a M\ri\Th + e\ 0) ) 
p 1 

-a M\rj\^ S in((l - p)a M\r]\^ + 0^)] 



Note that l^l 3 ^ ~ t is a decaying mode in the contract- 
ing phase and thus the last terms inside the square brack- 
ets in the above formulae can be neglected compared to 
the first ones. Since 



2(^°)^) , + d4°)4°) / ) 

l-c4 0)2 -M° )2 
-2rfVl 0) ' + M 0) 4 0) '), (70) 



then combining all these results we get: 
57 - Ci|t/|t=f . 



(71) 



For a\rj\ -Cm 1 ,M 1 , then differentiating (|55|) and 
56]) with respect to we obtain 



(0)' 1 - 3p ((1 - p)a nt) 2 |?7|i"p 



l-p 



a r(^) 



(72) 



l-3p 



6 (0)' _ l-3p((l-p)a AJ) 2 " 



i-p 



when p > 1 /3 and 



(74) 



when p < 1/3. Then we can solve Equation ([67]) to get: 

3-7p 1 

(57 ~ C 2 |7?| , p > - 

<5 7 - C 3 cos(fc|r/| +6»«5 7 ) . p<i (75) 

In the above expressions for £7, Ci, C2 and C3 are com- 
plicated prefactors in front of the ^-dependent terms. 

In summary, we see that the interactions give only a 
sublcading correction £7 to 70. 



D. Solutions for 



and 



Finally, let us consider the homogeneous component of 

(2) 

the second order fluctuations of the scalars, namely, cf)\ 

and 4> 2 ■ If we only consider the energy density up to 
second order, these second order field perturbations give 
a contribution through their coupling to the background 
fields. In the following we find the solutions of (|4"2l and 
study the effects of the induced terms in the stress-energy 
tensor on a possible bounce. 

Given the solution for the gauge field 7 obtained in the 
last subsection, it is easy to rewrite Eqs. (l4"2l as: 



,( 2 ) 



2™2„,(2) _ „(fcV- 7 ' 2 )„,(0) 



+ a m u 



,(2)" 2 Af 2„(2) 

t 9 + a ivi iio 



d (fcV-V a ) .,(o) 



(76) 



2 ' 
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where we made use of the fact that 

F^F^ = 2(fcV - a 2 f)/a 4 . (77) 

Equation (|76|) has the same form as the zero-th order 
equation but with a small source term generated by the 
interaction with the gauge field. This equation can be 
solved using the Born approximation (details are given 
in the Appendix). The general solution is the sum of the 
general solution of the homogeneous solution plus the 
solution including the source which has vanishing initial 
data. The inhomogeneous term is suppressed by the cou- 
pling constants c and d compared to the homogeneous so- 
lution, but, as shown in the Appendix, it scales as a high 
power of ?7 _1 . Via the coupling to the background scalar 
fields, the above second order terms enter into the expres- 
sion for the energy density to second order. The signs of 
the corresponding terms in the energy density are indef- 
inite in the sense that they depend on the phases of the 
initial field configurations. Since it is these terms that 
dominate the energy density near the bounce, we find 
that whether a bounce occurs or not depends sensitively 
on the phases in the initial conditions, and that in fact 
in the case of many plane wave modes initially excited, a 
bounce requires very special phase correlations. 

V. EVOLUTION OF THE COMPONENTS OF 
THE ENERGY DENSITY 

In the previous section we have solved all of the field 
equations up to second order in the amplitude of the fluc- 



tuations. We have found the scaling in time of each field 
at each order. Now we are ready to look at how all of 
the terms in the expression for the energy density p^°\ 
p' 1 ' and p^ at various orders in perturbation theory 
(namely, Eqs. p7|) - (fT9l ) scale in time. This analysis is 
straightforward but very important if we are to determine 
whether a bounce is possible, since in four space-time di- 
mensional classical Einstein Gravity with flat spatial sec- 
tions a bounce can only happen when the negative terms 
in the energy density catch up to the positive contribu- 
tions [n[. 

In the following we give a table of how each term con- 
tained in p scales with time as the background cosmol- 
ogy bounce point (the bounce which is achieved in the 
absence of radiation and scalar field inhomogeneities) is 
approached. We will identify the terms which dominate 
in this limit. This will give us a good indication under 
which conditions a bounce can occur. The tables are 
structured as follows: the first line "Terms", indicates 
which term we are considering, the next set of lines "Be- 
havior" gives the scaling in time in the various limits and 
in the two relevant ranges of the parameter p which in- 
dicates the equation of state, and the last line gives the 
sign with which the term contributes to the energy den- 
sity. Note that we focus on the growing mode solution 
to each field (which is constant for small r\ in the case 
p < 1/3). We give separate tables for terms of zero-th, 
first and second order in e. 

a) For terms contained in p^ : 



Terms 




m 2 < )2 








a 3 i (a\i]\ > mT 1 ) 


a~ 3 (a ?y ~> m -1 ) 


a 3 v (a\r\\ > m _1 ) 


a- 3 (a\r]\ > A/" 1 ) 


Behavior 


[ P>\ ] 


a -6+§C a l ? 7l < m_1 \ 
p>\ 


[ P>\ ' 


o _ 6+ |/oH 




p<\ 


n /a|?7| < mr 1 , 
[ P<\ ] 


p<h 


[ P<\ ' 


Sign 


Positive 


Positive 


Negative 


Negative 


Definite 


Definite 


Definite 


Definite 



T 



b) For terms contained in pW; 
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Terms 


X °U 1 
01 <P \ 


01 0i 


iOil 
-02 02 


Behavior 


a~ 3 ~p(\r]\ > Maxffc" 1 , (am)" 1 }) 


a_3 ( hi > Max{k~ 1 , (amy 1 }) 


a~ 3 ~t > Maxik' 1 , (aM)- 1 }) 


i i , r 7 — 1 / \ — li 

a -§~(H e t fc Slam) X K 


i i , r 7 — X ( \ — li 

( p>^ } 


I i , r 7 — 1 ( ~\/f\ — li 

a -|-5 (hi G t fc >( aM ) ]) 

P > 5 


1 i « i/< r T — 1 / \ — 1 1 

a -6( hi < Minjfc \ (am) ^ 
p > 3 


1 1 ^ 7i /T • r 7 — 1 / \ — 1 1 

a -6+f / hi < Mm{/s MH }\ 


1 1 ^ ji r - r j — 1 / li /r\ — 1 1 

-6/ hi <^Min{k , (aM) L }^ 
p > 3 




o _| /M e [fc- 1 ,(am)- 1 ]x 


Q Jr/le^.CaM)" 1 ]^ 


,|tj| < Mm{fc- 1 ,(am)- 1 }^ 


a o M <Min{k-\(am)- 1 }^ 


, M « Minfr 1 , (aA/)- 1 } v 


Sign 


Indefinite 


Indefinite 


Indefinite 



Terms 


-02 O) 0^ 




a" 3 (|r?| > Max{k-\ (aM)- 1 }) 




a - i+ xAri\G[k- 1 ,(aM)- 1 U 




P> s 


Behavior 


a -6+| / hi « Mm{fc-\ (aM)" 1 } v 
P > s 












a „ M^Mmik^^aM)- 1 }, 




P<k 


Sign 


Indefinite 



These terms, however, all vanish if the energy density is 
defined by spatial averaging. 



c) For terms contained in p( 2 '; 



J 



Terms 


0i^ 


01 01 




2 ,(1) 2 

m 0i 


m 2 0( o) 0i 2) 


Behavior 


a 3 i(k\r)\ > 1) 


a -7- ? (a 77 3> m _1 ) 


a- 5 (k\r]\ > 1) 


a- 3 (fc|r/| > 1) 


a _7+ p (oh| 3> m _1 ) 


o-6 ( *hl<<i) 


_ir / a ?? -C 771 -1 \ 
a ^ } 


o _ 8+f( %|<K 
P> S 




a -¥+i( a hl 

p > 5 


0( *N<i) 

P < s 


o / a hl < \ 
p < h 


p 3 


o ( *M<i) 

p 3 


a -¥+| ( a hl 

P < s 


Sign 


Positive 
Definite 


Indefinite 


Positive 
Definite 


Positive 
Definite 


Indefinite 



Terms 


-4^ 


-0 2 O) 0f 


^02 1)2 




-M^f^ 


Behavior 


a 3 * (fc|r/| > l) 


a 7 p (a|j7| > M _1 ) 


a- 5 (fc|r;| > l) 


a" 3 (fch| > l) 


a- 7+ p(a|?7| > M- 1 ) 


-«(*M<i) 


a 2( } 


a- 8+ t( fc|7?l< < 1 ) 
P> S 


a -e+f /fcN^l) 
P > 5 


o-¥+f( a hl« M_1 ) 


0( *N<i) 
p 3 


Q ,a\v\«M- 1 ^ 


a -2 (%!<<!) 
P < s 


o ( *N<i) 
p < 5 


a _ ¥ + , { a\r 1 \«M- 1 ) 


Sign 


Negative 
Definite 


Indefinite 


Negative 
Definite 


Negative 
Definite 


Indefinite 
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Terms 


„-4/„2„,2 i „-2-,2 

a k 7 +a j Q 


{-C(p\ -d<p 2 )(a k*y +a 7 ) 


„-4;„2„, r„, 


a 2 7(57 


Behavior 




a~ 7 (a\Ti\ > mr l ) 


aT 8+ p (a 77 > m -1 ) 


a~ 7 (a|?/| > m _1 ) 


a -io + i,a\r,\<m-K 




-10+2 / a t) < mT 1 \ 
P > 1 


4 /ffl|??| < m- 1 s 


a _ 4( a| ?7 |«m-^ 


4 /a|j7| < m- 1 s 


Sign 


Positive 
Definite 


Indefinite 
(Depending only on c and d) 


Indefinite 


Indefinite 



These terms to not vanish upon spatial averaging. 

Note that we have expressed the time dependence in 
terms of the dependence on the scale factor a(t). At 
this stage, we only need to focus on the exponent of the 
power-law scaling. The more negative the power is, the 
more rapidly the term grows in a contracting phase (since 
a(t) is decreasing with time). 

As mentioned earlier, the conditions for a bounce to 
occur in four space-time dimensional classical Einstein 
gravity with fiat spatial sections is that the total energy 
density reaches zero during the contracting phase. Thus, 
there needs to be a negative definite term which starts 
out small but grows faster than the positive definite terms 
due to the regular scalar field and regular radiation. In 
the absence of radiation and scalar field inhomogeneities, 
it is the contribution to the energy density of the ghost 
field <p2 which plays this role. 

From the table we see that there are three kinds of 
terms: positive definite, negative definite and indefinite 
ones. The first set contains the kinetic and potential 
terms of the normal scalar as well as the free energy 
density of the gauge field, the second set is made up of 
the kinetic and potential terms of the ghost scalar, while 
the third set contains terms which arise due to the cou- 
pling terms between scalars or between scalars and gauge 
fields. 

Looking first at the terms which are independent of the 
coupling term between the fields, we see from the first 
line of the "Behavior" set of lines in the third table that, 
indeed, in the presence of radiation the energy density in 
radiation grows faster than that in the two scalar fields, 
thus preventing a bounce. In the presence of coupling 
between the fields, however, there are terms which scale 
with a larger negative power of a(t). The signs of some 
of them, however, depend on the initial phases for the 
linear fields 7, <j>^ and tfi^ . 

Note that the signs of the scalar coupling terms are 
determined by the evolution of each field and thus are 
hard to be identified in a general analysis. The same is 
true for the gauge coupling terms (the last two in the 
third table). However, the coupling terms between the 
scalar fields and the gauge field (the third to last in the 
third table) can be made negative/positive definite easily 
by setting the signs of the coefficients c and d to be both 
positive/negative. 

It is reasonable to assume that the contracting phase 



begins with the regular scalar field dominating the energy 
density, and that the contribution of the Lee- Wick scalar 
is much smaller. For single Fourier mode initial condi- 
tions of the radiation field, this can be achieved with the 
appropriate choice of the initial phase (sec Example 1 in 
the following section containing our numerical results). 
However, for multiple initial radiation Fourier modes ex- 
cited any initial phase difference between the modes will 
produce a contribution with the wrong sign and will thus 
prevent a bounce (see Example 3 in the following sec- 
tion). In the presence of an infinite set of modes, the 
phase correlations required to obtain a bounce thus ap- 
pear to have negligible measure in initial condition space. 
Thus, even in the presence of coupling between scalar 
fields and radiation, the Lee- Wick bounce is unstable. 

The bounce, if it exists, will happen at a time which 
can be chosen to be t = 0. Its duration (the time interval 
lasting from the time the Hubble radius stops decreasing 
in the contracting phase until when it starts expanding 
in the post-bounce phase) will be denoted by At. Since 
the various components of the energy density scale with 
different powers of a(t), it is clear that the duration of 
the bounce will be shorter or equal to the Hubble radius 
Hmax (which gives the time scale on which the ratios 
of energy densities in different components change) at 
the beginning of the bounce phase. For the background 
bounce model, we have H m£iX ~ m. 

There are two kinds of bounces according to the dura- 
tion of the bounce phase i) If the period At ~ m _1 , the 
bounce will go from the time 



in- 



to the time 



t 



B + 



(At) 



(At) 



1 

2m 



1 

2m 



(78) 



(79) 



with a low speed. We call this a "slow bounce" . In this 
case, the universe will enter the bounce period at the 
critical time t c ~ m, and only the a\ij\ 3> m _1 approx- 
imate solutions of the previous tables will be applicable 
and not ones for the interval a, 1 77 1 <C m . ii) If the pe- 
riod At <C m , the bounce will happen in a very short 
time with very fast speed. This can be called the "fast 
bounce" . In this case, the universe evolves from the far 
past (— U with \ti\ ^> 1) to t = 0, passing through the 
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point t c ~ m, then entering into the region \t\ < m before 
finally reaching the bounce point. In this case, both of 
the two approximate solutions of the field evolution will 
be applied. 

Let us now consider the necessary conditions for a 
bounce (as we have indicated above and will sec from the 
numerical analysis, these conditions arc not sufficient - in 
addition to the conditions which follow, appropriate cor- 
relations in the initial phases are required). We start in 
the region of time a\r]\ 3> mT l . We study the conditions 
required to have the terms that might give a bounce grow 
relative to the other terms during this phase. If the con- 
ditions are not satisfied, or the bounce does not happen 
even if the conditions are satisfied, then a bounce may 
still occure in the a. 1 77 1 <C to -1 region. The conditions 
for the terms in the energy which could compensate the 
positive radiation contribution to become dominant are 
then studied. If these conditions are not satisfied, either, 
then a bounce is impossible. 

A necessary condition for a bounce to be possible re- 
quires the growth rate of one of the indefinite sign terms 
in the third table above exceed all that of all of the pos- 
itive definite terms. In the a\rj\ 3> mT 1 region, this re- 
quires 8 — l/p > 5, which equivalent ly constrains the 
equation of state parameter w to be in the range w < 1. 
If this condition is satisfied in this region, then a slow 
bounce may happen depending on the choice of the ini- 
tial phases. 

If the condition is not satisfied in the a\rj\ 3> to -1 
region, the universe may evolve into the o.| 77 1 -C mT 1 
region, in which the evolution of the fields are different, 
and new constraints on p and w will arise if a bounce is 
to be possible. Following the above logic, we find that 
the conditions under which a bounce might happen are 
much looser, namely w > —7/6. 

To summarize this section: we have identified neces- 
sary conditions for a bounce to occur. Whether one 
actually does occur even if the conditions are satisfied 
depends on the initial phases of the fields. This must 
be studied numerically. In the following section we will 
give one example of specially chosen phases for which a 
bounce is possible. However, when we look at a more 
general choice of phases, the bounce will not occur. 



VI. NUMERICAL RESULTS 

In order to support the analysis in the last section, we 
performed numerical calculations. Such numerical work 
is necessary because our analytical analysis is only ap- 
proximate. In particular, we worked in perturbation the- 
ory up to order second order in e. In addition, even in 
cases where our analytical analysis would indicate the 
possibility of a bounce, the pcrturbative analysis will 
break down near the bounce point, and there is no as- 
surance that the trends seen in the pcrturbative analysis 
will persist. 

We have numerically solved the full nonlinear equa- 



tions of motion for the matter fields in the presence of a 
homogeneous expanding background cosmology. The ho- 
mogeneous cosmology is obtained numerically by solving 
the first Friedmann equation 



where G is Newton's gravitational constant (related to 
the Planck mass used earlier), and p is the total energy 
density, averaged over space. 

Figures [T][5] are two groups of numerical results with 
different parameters. In both cases we choose the initial 
energy density of the gauge fields to be larger than that 
of the Lee- Wick scalar, but less than that of the normal 
scalar. These initial conditions correspond to the situ- 
ation we are interested in, namely starting in a matter- 
dominated contracting phase in the presence of some ra- 
diation which is sub-dominant. Fig. [TJ [2] and [3] show an 
example with parameters c > and d > 0. We choose 
initial conditions in which a single Fourier mode fluctu- 
ation is excited, and in which the phases are chosen as 
indicated in the figure caption. For these initial phases, 
we obtain a bounce. In Fig. [TJ we see that the equation 
of state w begins with a value slightly larger than 0, and 
then evolves to some nearly fixed value. For the case of 
our initial condition choice, it appears to be w ~ —0.6, 
in the region where the bounce is allowed to happen. 
At the bounce point, the equation of state will drop to 
—00, while after the bounce, the equation of state will 
rise again to w ~ 0.6. Fig. [2] is the plot of the scale 
factor in this case which shows explicitly the occurrence 
of the bounce. Fig. [3] gives a comparison of the en- 
ergy densities of some components during the process. 
Initially, we set the energy density of the gauge field 7 
to be between the normal scalar and Lee- Wick scalar. 
When the evolution of the universe enters into a region 
with nearly constant w, the gauge-coupling component 
of energy density p c will grow very fast. It is negative 
and thus enables the negative part of the energy den- 
sity to catch up with the positive one, thus allowing the 
bounce to happen. For the inhomogeneous fluctuation, 
we choose the wavenumber to be k ~ O.OlhMpc^ 1 which 
corresponds to a scale which is observable by CMB and 
LSS experiments. 

Figs. |H [5] and |6] give the corresponding results in the 
case when we choose c > while d < (with all initial 
conditions identical). This case seems dangerous because 
the contribution of the Lee- Wick scalar to the fluctuation 
terms could lead to an instability. However, as we have 
mentioned before, since the effects of the Lee- Wick scalar 
are less than that of the normal scalar, it is still possible 
for the bounce to happen. Fig. [T] shows the equation of 
state of the system. We can see that the evolution of w is 
about the same as that in Fig. [TJ since the change of the 
sign of d does not alter the result too much. Fig. [S]is the 
behavior of scale factor in this case while Fig. [BJgives the 
comparison of the energy densities of all components. 

A change in the phase of the initial radiation field ve- 
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FIG. 1: The evolution of the equation of state w w.r.t. cosmic 
time t (horizontal axis), in the first simulation, a simulation 
with only a single Fourier mode excited and phases chosen 
as indicated below. We see that w drops to — oo, indicating 
that there is a bounce. The background fields are plotted 
in dimensionless units by normalizing by the mass m rec — 
W~ 6 mpi while the time axis is displayed in units of m^} c . 
The mass parameters m and M were chosen to be m = 5m rec 
and M = 10m rec . The initial conditions were chosen to be 
ji ~ -1.85 x 10 5 m rec , ji ~ 7.35 x 10 6 m 2 ec , tj> u ~ 1.015 x 
10 5 m rec , 0ii ~ 6.39 x 10 5 m 2 ec , (j> 2i ~ 2.54 x 10 2 m rec , fai ^ 
—4.96 x 10 3 m 2 . ec . The coefficients c and d are chosen to be 
c = 10~ 10 M 2 ec and d = 10~ 10 'M r 2 K . The wavenumber is 
fc ~ 0.01/iMpc" 1 . 




-1.0 -0.5 0.5 



FIG. 2: The scale factor of the universe in the same simulation 
that leads to the evolution of the equation of state shown in 
Fig. [T] From the plot we see that the bounce happens at 
t = 0. 



locity will not change the results (if we keep the other 
initial conditions fixed). On the other hand, if we flip 
the sign of the initial velocity of one of the two scalar 
fields, then the sign of the dominant contribution to the 
energy density as we approach the bounce will flip and 
this will prevent a bounce. If we use initial conditions 
containing two excited Fourier modes, then we obtain a 
bounce only if the signs of the initial field velocities are 
both chosen as in the first run whose results are shown 
here. Different phases for the scalar field velocities of the 
two modes destroys the possibility of obtaining a bounce. 




100000 I — I ' 1 ' 1 ' L_ 

-1.0 -0.5 0.5 

t 



FIG. 3: (Color online.) Energy densities of </>i, 02 and 7 in 
the system with parameters chosen as in Fig. [T] The curves 
from up to down are: p ( p 1 (black), p^ 2 (red), p c (dark cyan) 
and p 1 (blue), respectively. The variables are also normalized 
with the mass scale m rec = 10 -6 mp(. 
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FIG. 4: The evolution of the equation of state w as a function 
of cosmic time t (horizontal axis). The behavior that w drops 
to —00 indicating that a bounce takes place. The background 
fields are plotted in dimensionless units by normalizing by 
the mass m rec — 10 -6 mp; while the time axis is displayed in 
units of m~l r . The mass parameters m and M were chosen to 
be m = 5m rec and M = 10m rec . The initial conditions were 
chosen to be ji ~ — 1.85x 10 5 m rec , 7; ~ 7.35 x 10 6 m^ ec , <f>u — 
1.015 x 10 5 m rec , ~ 6.39 x 10 5 m 2 rec , cf> 2i ~ 2.54 x 10 2 m rec , 
4>2i — —4.96 x 10 3 m 2 ec . The coefficients c and d are chosen 
to be c = fO _10 M 2 ec and d = -f 0~ 10 M 2 ec . The wavenumber 
k ~ 0.01/iMpcT 1 . 



Figures [S] and O show the results for the equation 
of state parameter w, the Hubble parameter H and the 
contribution of the various components to the total p 
in the case of a simulation in which two Fourier modes 
are excited, with velocities of both scalar fields having 
opposite signs from those in the previous example. As is 
obvious, a Big Crunch singularity occurs. 
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FIG. 5: The scale factor of the universe driven by the system 
with parameters chosen as in Fig. [4] From the plot we see 
that the bounce happens at t — 0. 
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FIG. 7: The evolution of the equation of state id as a function 
of cosmic time t (horizontal axis). The background fields are 
plotted in dimensionless units by normalizing by the mass 
m rec — 10~ 6 mp; while the time axis is displayed in units 
of ni^ec- The mass parameters m and M were chosen to be 
m = 5m rec and M = 10m rec . This plot is the evolution 
of the system with two Fourier modes combined together. 
The one is of which the wavenumber k ~ 0.01hMpc~ with 
initial conditions 7, ~ —1.85 x 10 5 m rec , 7, ~ 7.35 x 10 6 m^ ec , 
<j>u ~ 1.015 x 10 5 m rec , 0ii ~ 6.39 x I0 5 m 2 rec , foi ^ 2.54 x 
10 2 m rec , 4>2i — —4.96 x 10 3 m^ ec , which if taken alone will 
give the bounce as has been shown in the previous example. 
The other is of which the wavenumber k ~ Q.04hMpc~ with 
initial conditions of the same initial values of the fields but 
the opposite signs of the scalar field velocity. From this plot 
we can see that the combination of the two Fourier mode 
will (generally) cause w blow up, thus preventing the bounce. 
This means that the bounce requires special fine-tuning of the 
initial phases for each Fourier mode. The coefficients c and d 
are chosen to be c = 10" 10 M£, C and d = -l(T 10 M£, o . 



FIG. 6: (Color online.) Energy densities of </>i, 02 and 7 in 
the system with parameters chosen as in Fig. 2] The curves 
from up to down are: p ( f >1 (black), p^ 2 (red), p c (dark cyan) 
and p 7 (blue), respectively. The variables are also normalized 
with the mass scale m rec = 10~ 6 mpi. 



VII. CONCLUSIONS AND DISCUSSION 

In this paper we analyzed in detail the possibility of 
obtaining a cosmological bounce in a model which corre- 
sponds to the scalar field sector of the Lee- Wick theory 
coupled to relativistic radiation. It is known that the 
scalar field sector of the Lee- Wick theory in the absence 
of other fields can yield a cosmological bounce [l[ . In fact, 
the universe will scale as non-relativistic matter with 
< w >~ both before and after the bounce. Thus, this 
model is a possible realization of the "matter bounce" 
scenario. However, this background is unstable to the 
introduction of radiation since in the contracting phase 
the growth of energy density in radiation will exceed that 
of matter and will lead to a Big Crunch singularity As 
has been shown in previous work [i~5j . the introduction 
of a Lee- Wick partner to radiation does not prevent this 
instability. In this paper, we introduced an interaction 
between the radiation field and the scalar fields. The 
interaction could help drain energy from the radiation 
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FIG. 8: The Hubble constant of the universe driven by the 
system with parameters chosen as in Fig. [7] From the plot 
we see that there is a singularity at t = 0. 



field to the Lee- Wick scalar, and thus could prevent the 
radiation from growing too fast to destroy the bounce. 

We analyzed the equations describing the evolution of 
the three matter fields (regular scalar field, its Lee- Wick 
partner and the radiation field) on a cosmological back- 
ground both analytically and numerically. Our analyti- 
cal analysis was perturbative and made use of the second 
order using Born approximation. The expansion param- 
eter is set by the initial amplitude of the gauge field. We 
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FIG. 9: (Color online.) Energy densities of (pi, 02 and 7 in 
the system with parameters chosen as in Fig. [7] The curves 
from up to down are: p ( f, 1 (black), p^ 2 (red), p c (dark cyan) 
and p~i (blue), respectively. The variables are also normalized 
with the mass scale m rec = 10~ 6 mpi. 



Appendix A: Green's function determination of u\ 

The solution for the second order scalar field correction 

(2) 

u\ can be determined using the Green function method. 
The general solution of (|76p is the sum of the solution 
Mo (77) of the homogeneous equation which solves the same 
initial conditions as u and the particular solution 8u(jf) 
wich vanishes at time r\i . The particular solution is given 
by 

rv 

Su(rj) = u a {r]) / drj ' e{r)')u b (r)')s(r)') 

J r)! 

-u b {ri) [ dr,'e{r)')u a (ri') S (ri') , (Al) 

where u± and ui are two independent solutions of the 
homogeneous equation, e(jj) is the Wronskian: 



solved the equations of motion for each field at each or- 
der, and obtained their approximations in different cases. 
We compared their contributions to the total energy den- 
sity, and derived necessary conditions for a bounce to 
happen. To support our analysis, we also performed nu- 
merical calculations. 

Specifically, we investigated initial conditions in which 
one or two Fourier modes of the radiation field and the 
scalar field fluctuations are excited. We found special 
initial conditions which indeed lead to a non-singular 
bounce. Changing the sign of the initial scalar field ve- 
locity will destroy the bounce solution. In the presence 
of two Fourier modes, we found that a bounce requires 
identical initial phases for the two modes. For general ini- 
tial conditions, we conjecture that the measure of such 
initial conditions which lead to a bounce is very small. 
We thus find that the addition of coupling terms between 
the scalar fields and radiation cannot save the Lee- Wick 
bounce background from the instability problem with re- 
spect to the addition of radiation (nor, for that matter, 
with respect to scalar field fluctuations). The instability 
problem with respect to anisotropic stress will be even 
worse. 



e(»7) = (u a u b -u b u a ) , (A2) 

and s(r]) is the source inhomogencity. 

Recall from the main text that the second order field 
correction terms satisfy the equations: 

(A3) 

(2) 

We will demonstrate the analysis for the case of u\ ' . 
Let us consider evolution for a short interval of time start- 
ing at some initial time r/i. Then, we can neglect the 
expansion of the universe in the equation of motion and 
take a(rf) = a(r)j). We are then interested in how the 
result scales in r\i. Using this trick, the solutions of the 
homogeneous equation can be taken to be 

u a (r)) = cos(w ln ?7) 

u b {rj) = sm(uj rn r]) (A4) 

and the Wronskian is 



e {v) = where u m = V a 2 m 2 . (A5) 
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Using the result for the background 7 from the main 
text, the source term becomes 



s~<{ri) = c 



(fcV-y 2 ) (o) 



D ~y\'U — *-* 9 ^1 

or 



ar 2 k 2 \r]\^^ x (A6) 
cos(2fc|ry| + 20 7 ) cos((l - p)am\r}\ + ^ 0) ) , 



since 70 ~ cos(fc|?7| + 6* 7 ). 
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Combining these results we obtain Note that if we only care about their scalings with re- 
spect to conformal time rj or scale factor a(t), the above 

(2) / \ f v d-qk i v solutions can be reduced to: 

u\ cos(w m ?7) / |?7|^tt x 

Jrji ^ \t)Tfl 

sin(w m 77)cos(2fc|77| + 26' 7 )cos((l-p)am|77|+6'f ) ) 9p _ 2 _ a 

fV dllk 2 p wijoc \q\*G=v oca 2 + p , (A9) 



+ sin(w m 77) / 1^315=17 x (A7) 

cos(w m ?7) cos(2/c|77| + 20 7 ) sin((l — p) am 1 77 1 + 9^) , 



and for the case of a matter dominated era where p = 2/3, 

4 2) ~ - costoJMV) f 1 -$^-H^ x lt is straightforward to show that 

J m a 6 {t)M 

sin(w M ?7) cos(2fc|7?| + 20 7 ) cos((l - p)oM|r/| + 0^ o) ) 

. , , r d v P . ... cx H- 6 tx a- 3 . (A10) 

J VI a 3 (t)M 

cos(lu m i]) cos(2fc|?y| + 20 7 ) sin((l - p)aM\rj\ + 6 { ° ] ) . 
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